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Abstract. Given two probability measures fj, and v we consider a mass trans- 
portation mapping T satisfying 1) T sends fi to u, 2) T has the form T = 
ip | ^^ | , where ip is a function with convex sublevel sets. We prove a change of 
variables formula for T. We also establish Sobolev estimates for ip, and a new 
form of the parabolic maximum principle. In addition, we discuss relations to 
the Monge-Kantorovich problem, curvature flows theory, and parabolic non- 
linear PDE's. 
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1. Introduction 

In this paper we study a class of mass transportation mappings having the form 

T = tp W\ 

with some potential f. The mappings of this type have been introduced in [5], [5]. 
Assume we are given a couple of probability measures fi — p dx and u = p\ dx. 
It has been shown that, under general assumptions, there exists a unique ip with 
convex sublevel sets A t = {x: ip(x) < t} such that 

T: x — > f(x) ■ n(x), 

where n(x) is the normal vector to dA t at x with t — <f{x) and T satisfies the 
equality v = /.io T^ 1 . We point out that the restriction of T to every level set dA t 
coincides (up to the factor t) with the Gauss map of dA t . In what follows we use 
the name "Gauss mass transport" for T. 

Mappings of this kind are closely related to several areas of research. They can 
be considered as "parabolic" analogs of optimal transportation mappings, which 
attract attention of researchers from the most diverse fields, including probability, 
partial differential equations, geometry, and infinite-dimensional analysis (see [36], 
[37], and [7]). In addition, they arise naturally in the Gauss curvature flow theory. 
Concerning transformations of measures of other related types, see [S], [5], |10) . 

The main goal of this paper is to establish some regularity properties of the 
mapping T. More precisely, we prove that T satisfies a change of variables formula, 
which can be considered as the weakest regularity property of T . 

The corresponding result in the elliptic case (optimal mappings between mea- 
sures with densities always satisfy a change of variables formula) belongs to McCann 
[29] , This result turns out to be quite useful for different applications. Applications 
of the change of variables formula include, for instance, the so-called above-tangent 



formalism which is a crucial technique in variational problems, PDE's, and proba- 
bility (see |37], 0, H). 

The paper is organized as follows. 

In Section 2, we briefly describe the main results of [9] that are used throughout. 
These are the results on existence and uniqueness of Gauss maps, a description of an 
important scaling procedure, and certain duality relations. In addition, we describe 
the relations to curvature flows and the parabolic Monge-Ampere equation. 

Our main result is proved in Section 3. We show that T satisfies the following 
change of variables formula: 

pi°T-J = p with J = det D a T, 

where D a T can be understood as the absolutely continuous part of the distributional 
derivative of T. One has 

J = i P d - 1 \D a ip\K, 

where |-D a </j| is the absolutely continuous part of the full variation of the vector- 
valued measure and K is the Gauss curvature of dA v i x y 

In Section 4 we establish some natural Sobolev a-priori estimates for (p. We 
emphasize that ip is not Sobolev but only BV in general. Under assumption that 
Pi = r d-i we show that for every p > 

C p ,r [ |V^| P+1 dp < f ^ P+1 dp+ f K-P P P Q +1 dU d -\ 

J A J A PO JdA 

Another natural question arising in the study of the Gauss mass transport is the 
validity of some parabolic analogs of the maximum principle. Applying the mass 
transportation arguments one can establish (see Section 5) the following form of the 
parabolic maximum principle: every smooth function / on a convex set A satisfies 
the inequality 



sup/ < sup/ + C(d) f \Vf\Kdx, 

A OA JC-f,i 



where 

C- f ,i ={x:xeA t nd conv(At)}, A t = {-/ < t] 

is the set of contact points for the sublevel sets of — /, conv(A t ) is the convex enve- 
lope of A t , and K is the corresponding Gauss curvature. This estimate is naturally 
related to the Gauss mass transport and the second-order nonlinear parabolic dif- 
ferential operator / i-> \V f\K (similarly to the Monge-Ampere operator in the 
classical maximum principle). The inverse mapping S — T _1 is associated with 
another parabolic differential operator: 

/ r .dct(/-Id + D e 2 /) 

where Dgf is the Hessian on S^ -1 . The corresponding maximum principle is proved. 

In Section 7, we are concerned with the regularity of the parabolic Monge- 
Ampere equation. In particular, we briefly explain how the arguments employed in 
[34) can be extended to our situation to prove Holder's regularity of ip. Thus we 
establish Holder's continuity of if assuming that pi,p2 € C 2 ' a (A) and dA is smooth 
and uniformly convex. 

The author express his gratitude to Vladimir Bogachev for valuable suggestions 
and remarks. 
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2. Existence and basic properties 



In what follows we denote by T-L m the m-dimensional Hausdorff measure on M. d , 
m < d. For Lebesgue measure we also use common notation A. We denote by S^ -1 
the unit sphere in R d (and by S^T 1 its upper- half). We also use the symbols Dg, 
Dg for the gradient and the Hessian on S 1 " -1 . 
It will be assumed throughout the paper that 
Al) the measure \i is supported on a compact convex set A 
A2) the measure v is supported on Br = {x: \x\ < R} for some R > 
A3) the measure fi is absolutely continuous with respect to A|,4 and v is abso- 
lutely continuous with respect to A|b r - 

We start with a brief outline of two areas of research closely related to the Gauss 
mass transport. 

1) Optimal transportation. 

Optimal transportation can be described as a problem of optimization of a certain 
functional associated with a pair of measures. The quadratic transportation cost 
W^in, v) between two probability measures /i, v on WL d is defined as the minimum 
of the Kantorovich functional: 



(1) m i-> / |xi— X2I dm(x\,X2), mcP(/i,y), 

where V(^jl, v) is the set of all probability measures on M. d x K d with the marginals /1 
and v\ here \v\ denotes the Euclidean norm of v € M. d . The problem of minimizing 
is called the mass transportation problem. In many cases there exists a mapping 
T: M. d — > R d , called the optimal transport between /i and v, such that v = /10T" 1 
and 

W2(n,v)= I \x-T(x)\*n(dx). 



If and v are absolutely continuous, then, as shown by Brenier and McCann (see 
[36]), there exists an optimal transportation T which takes fi to v. Moreover, this 
mapping is /i-unique and has the form T — VVF, where W is convex. Assuming 
smoothness of W, one can easily verify that W solves the following nonlinear PDE 
(the Monge- Ampere equation): 



pi(VW) 

In fact, this equation is satisfied in a certain sense without any smoothness assump- 
tions (see Section 3). 
2) Geometric flows. 

We refer to [TB], [T7] for an account in geometric flows. Let {r f } C M. d be a 
family of embedded hypersurfaces. Denote by V{x, t) the velocity in the direction 
of the inward normal — n(x) at a point x £ T t . We say that {T t } satisfies a surface 
evolution equation (or {r f } is a geometric flow) if V satisfies 

(2) V = f(x,t,n,-Dn) 

for some given function /. If / — H is the mean curvature, then T t is called the 
mean curvature flow. If / = K is the Gauss curvature, then T t is called the Gauss 
curvature flow. 

The Gauss curvature flows have been introduced by Firey [T5] as a model of 
the wearing stone on a beach. The existence and uniqueness of a Gauss curvature 
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flow starting from a smooth initial convex surface has been obtained by Tso [34) by 
solving a corresponding parabolic Monge-Ampere equation. He proved, in partic- 
ular, that r t remains convex and shrinks to a point in finite time. The same result 
for the mean curvature flow has been obtained by Huisken [21]. More on Gauss 
curvature flows see in [3]. 

The main problem arising in respect with non-convex initial surfaces is the even- 
tual singularity of the solution. It turns out that in general T t becomes singular in 
finite time. To overcome this problem several notions of generalized solutions have 
been proposed. A weak notion of a solution to (J2J) have been introduced by Brakke 
[11) . He proved the existence of the mean curvature flow for any initial data in 
some generalized measure-theoretical sense. According to the level-set method (see 
[T7]), the family {F t } is considered as level sets of some function u(t, x) satisfying a 
nonlinear parabolic equation in viscosity sense. Finally, it is known that sometimes 
the solutions to curvature flows can be obtained as scaling limits of certain ellip- 
tic or parabolic equations. For instance, the mean curvature flow can be obtained 
as a singular limit of the solutions to Allen-Cahn or Ginzburg-Landau equations 
(see [22], [31]). It has been shown in [9] that Gauss curvature flows starting from 
convex surfaces are singular limits of some optimal transportation problems. More 
precisely, the following result has been proved in [S]. 

Theorem. Let A C M. d be a compact convex set and let p = po dx be a 
probability measure on A equivalent to the restriction of Lebesgue measure. Let v — 
Pi dx be a probability measure on Br = {x: \x\ < R} equivalent to the restriction of 
Lebesgue measure. Then, there exist a Borel mapping T : A — > Br and a continuous 
function ip: A — > [0,R] with convex sub-level sets A s = {ip < s} such that v = 
fi o T _1 and 

T = ip ■ n H d -almost everywhere, 

where n = n(x) is a unit outer normal vector to the level set {y: y>(y) = <p(x)} at 
the point x. 

If ip is smooth, the level sets of ip are moving according to the following Gauss 
curvature flow equation: 

(3) x(s) = -s d - 1 P±M K (x)-n(x) 

Po(x) 

where x(s) € dAn- s , < s < r, x(0) € OA is any initial point satisfying ip(x(0)) = 
R. 

Remark 2.1. 1) The theorem does not guarantee that the boundary dA is 
exactly the level set {ip — R}. Nevertheless, one can easily check that this 
is indeed the case when A is strictly convex. 

2) It is not clear in general whether {x: ip(x) — 0} contains a unique point or 
just has Lebesgue measure zero. 

3) The case pi = po — C corresponds to the standard Gauss curvature 
flow. The asymptotic behavior of dA r for small values of r is a standard 
problem in differential geometry. For the classical Gauss flow it is known 
that dA r is asymptotically spherical in shape for values of r close to (see 
[3]). This problem has not been studied so far for the flows of the type ([3]). 

4) Potential (p is not Sobolev in general, but admits a bounded variation (BV) 
(see [1]). The distributional derivative of ip can have a singular component 
in the n-direction. 
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In addition, T is unique and admits an inverse T _1 (see [5], Section 3). 

Let us briefly describe the idea of the proof and some important related facts. 
The potential p is a pointwise limit of a sequence of functions {ift} with convex 
sublevel sets. To construct ift we consider the optimal transportation VWf of \x to 
v o , where St(x) = x\x\ l . Let us set 

|VWt|TTt 

Clearly, Tt pushes forward /i to v. Choose Wt in such a way that mh\ x ^A Wt(x) = 0. 
Dehne a new potential function ift by 

1 



Then one has 



It was shown in 9 that 



T t = <p t - 



lim ipt n — ip, lim T tri = T 

almost everywhere (for a suitable subsequence i„ — > oo). 
The dual potentials 

W7(y)= sup((x,y)-W r t ( a ;)) 

of the corresponding dual Monge-Kantorowich problem define via renormalization 
another natural convergent sequence 

It was shown in [9] that 
pointwise, where 

H(r, 6): B R = [0, fl] x S^ 1 -J- R, 
H(r,9) = sup (0,sg). 

We warn the reader that in [5] we deal with a slightly different potential iftt — 
H t (r,6)r. 

Let us describe the expression for T t ~ 1 in terms of H t . To this end we fix n e S d ~ 1 
and introduce local coordinates • • • , on in a neighborhood of n. We 

assume everywhere below that 

dn 

constitute an orthonormal basis in the tangent space of S^ 1 at n. Then the 
following relation holds 

d-i 

TT\y) = (H t + J^(Ht)r) ■ n + £(#tk • e*. 



In the limit t — > oo one has 

d-l 

T-\y) = H ■ n + J^-ff* • e< = H ■ n+ D e H. 

i=l 

Remark 2.2. In what follows we often choose the following convenient local coor- 
dinate system on .S d_1 . We take the center of «S d-1 for the origin and introduce 
the standard Euqlidcan coordinates in M. d such that n becomes the North Pole: 
n = (0, 0, . . . , 1). A neighborhood of n can be parametrized by 



(01, . . . , 9d-i) — > ■ ■ ■ , #d-ij (1 — 
In particular, one has at n: 



zV /2 



Clearly, (r, 6\, . . . , 9d-i) is a parametrization of a cone with the vertex at the origin. 

Now we describe the relation between the Gauss mass transport and the para- 
bolic Monge- Ampere equation. 

Several parabolic analogs of the elliptic Monge-Ampere equation have been in- 
troduced by Krylov (see |25j). He also proved some forms of the parabolic maximum 
principle (see also [35]). 

Let [i = po dx be a probability measure on an strictly convex set A. Consider a 
Gauss mass transportation 

= Wp~\ 

sending fi to a measure v = p\ dx on Br : = {x : \x\ < R}. 

Example 2.3. Assume d — 2 and fix a standard coordinate system [x\, x%). Assume 
that the functions below are smooth. Introduce the global polar coordinate system 
(r, 6). One has 

T- 1 = H-n + H e -v, 

n = (cos 0, sin 6 1 ), v = (— sin^, cos^). 

Let us compute the derivative of T -1 in polar coordinates: 

I7 1 = H r ■ n + Her ■ v 

Tg 1 = H -n e + H e ■ n + He ■ v fl + H ee • v = (H + H ee ) ■ v. 

Taking into account that detD(r,0) = ± one gets dctLT" 1 = ME±EsA _ Finally, 
by the change of variables formula 

(a\ (tj i v \ H r{H + Hee) 

(4) pi = p (H ■ n + H g ■ v) . 

r 

Let us describe a standard trick which allows to rewrite ^ in the form of the 
parabolic Monge-Ampere equation. Introduce another variable on X2 < 0: 

Z = -Ctg 9, 7T < < 2tt. 

Thus = arcctg(— z). Instead of iJ it is convenient to work with 

u = \[\ + z 2 H. 
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Note that u is just the restriction of the corresponding 1-homogeneous support 
function Ha t with a fixed r to the line X2 = — 1. In particular, u is convex in z. 
Taking into account that Jj- = j^rp ~§g > one can easily compute 



zH + Hg H + H ee 

u z = — , u z 



VT+Z* (1 + Z 2 )5 

Finally, we set 

T = T" 1 o (r,arcctg(-z)). 
Writing this mapping in coordinates (xj^JCa) as a function of (z,r), one gets 

T -1 = (u z ,zu z -u) = (u z ,u*(u z )), 
where u* is convex conjugated to u with respect to z-variable 



drdz. 



u*(z,r)— sup (xz — u(z, r)J . 
lei 1 

The change of variables formula takes the form 

5 Ur . Uzz = — r— ^pi (-7=5 , («,r et + xl. 

Note that (|5|) can be considered as a parabolic Monge- Ampere equation. 
In addition, ([5]) can be easily interpreted from the point of view of mass trans- 
portation. Indeed, let us set 

r ( rz, —r 

Then v is a measure on R x [0,R] which coincides with the image of v under the 
mapping 

(rz, -r) 
(x,y) i — > ■ 

Further, /i is the image of v under T _1 . Function u is convex in z and increasing 
with respect to r. 

All these computations can be generalized to the multidimensional case. One 
has T- 1 = H • n + Y^=l H er e t = H • n + DgH and 

„ ,^ 11, -dct(H- Id + D*H) 
(6) p 1 =pa(T ) . 

Here DgH denotes the Hessian of H on the unit sphere. For computing DgH it is 
convenient to deal with the local polar coordinate system as described above. In 
this case DgH(n) can be represented just by the matrix (dg .g .H). Note that 

det(H ■ Id + D 2 e H) = * H r 1 



K{T- r y ' \V<p(T- l )\' 
Finally, let us define coordinates (z,r) and the corresponding chart 

V{z,r): {-R < x d < 0} ->■ B R , 
(xx,...,Xd) = (zi,...,z rf -i,-l) = V{z,r). 

V 1 + 4 + ■ ■ ■ + 4-1 

Now we introduce a new potential u 



u= x /l + z 2 -\ Vz\_ x S 
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and verify the following proposition by direct computations. 

Proposition 2.4. Assume that T is smooth. The following representations hold 
on —R < Xd < 0: 

1) 

T" 1 = (V,«, (z, V z ii) - u) = (W z u,u*(W z u)), 

where 

u*(z,r) — sup z) — u(z, r)), 



2) 



T" 1 = T- 1 o K(z,r). 



det(H -ld + D 2 e H) = (i + z 2 + ... +z a_ i j|(d-i) det£ ^ 



3) ifte change of variables takes the form 

l2„. Pi 



• det £>fu = 



where 
Pi = 



l + z 



'*d-\ 



po(V z u, (z, V z u) - u) ' 

(- 



d-l 



rzi,-- ■ ,rz d -i,-r 



l + zl + 



More on the parabolic Monge-Ampere equation see in Section 7. 



3. Change of variables 

Let A be any convex compact set of positive volume and let T: A — > B be a 
Gauss mass transport between two given probability measures and v satisfying the 
assumptions specified in the introduction. To prove the change of variables formula 
for the Gauss mass transport we need to define the Gauss curvature for sufficiently 
"large" amount of points x S dA. To this end we consider the corresponding 
support function 

H A {6) = suv{6,x). 

Here we assume that 9 € M. d . Clearly, Ha is Thomogeneous and convex. Hence, 
by the Alexandrov theorem Ha is almost everywhere twice differentiable. Recall 
that every convex function V is a.e. twice differentiable in the Alexandrov sense, 
i.e. for almost all x there exists a matrix D^V(x) (the absolutely continuous part 
of the second distributional derivative) such that 

(7) \V(y)-V(x)-(VV(x),y-x)-±(D 2 a V(x)y-x,y-x)\=o(\y-x\ 2 ), y^x 
(see [H]). 

Remark 3.1. A parabolic analog of the Alexandrov theorem for monotone-convex 
functions was proved by Krylov (see }25j). 

Definition 3.2. In what follows we say that /: M —¥ M, where M is a Borel set is 
differentiable at x £ M in the sense of Alexandrov if (J7|) holds for y G M, 
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In particular, homogeneity implies that for every fixed r > the function HA\dB r 
is twice differentiable for H d_1 -almost all x £ dB r . 
Recall that H is defined as follows: 

H(r,6) = sup (9,x). 

ees d ,xeA r 

Lemma 3.3. For ^-almost all x e A and all < r < R the function H\qb t is 
twice differentiable at r ■ n(x) in the Alexandrov sense. 

Proof. It was noted above that H\oB r is twice differentiable in the Alexandrov sense 
for H d ~ 1 -almost all y € 8B r . Hence by Fubuni's theorem the set of all y such that 
H\dB r {y), r — \y\ is not twice differentiable in the Alexandrov sense has ^-measure 
zero. The claim follows from the fact that T pushes forward fi to v. □ 

Next we want to define the Gauss curvature for an arbitrary convex surface dA 
% d_1 -almost everywhere. Let us recall how the Gauss curvature can be defined in 
the smooth case. 

Assume that dA is a level set of some smooth function F. Then 

VF(x) 

n(x): = Wf&T\' 

Let {ei, e2, • ■ • , e<j-i} be an orthonormal basis such that n_Lej for every 1 < i < d— 1. 
Then 

K{x) = dctDn{x), 

where D is the differential operator on the tangent space to dA. Computing this 
expression in Euclidean coordinates one gets 

K — det Aij , where 

(8) A^(x) = 1 (D 2 F(x)e i ,e j ), 1 < i,j < d- 1. 

In particular, this formula is applicable when the surface is represented locally as 
the graph of a convex function 

F = W(xi, . . .,Xd-i) - x d , 

where {e,} can be obtained by an orthogonalization procedure from a basis tangent 
to F at some point. 

It is convenient to compute the Gauss curvature in terms of the support function. 
The following lemma (well known for smooth surfaces) gives another practical way 
of computing. 

Lemma 3.4. Let dA be a convex surface which coincides with a graph of some 
convex function W(xi, . . . , x^-i) in a neighborhood ft of x and n(x ) is unique 
at xq . Then the following are equivalent 

1) W is differentiable at x in the Alexandrov sense and D^W(xo) is nonde- 
generate, 

2) H\s d _ 1 is differentiable at n(xo) in the Alexandrov sense and 



det (if • Id + (D 2 e ) a H^j o n (a;o) ^ 0. 



Proof. Choosing an appropriate coordinate system, we may assume without loss of 
generality that xq — 0, W(0) = and VW(0) = 0. Let us assume for a while that 
W is smooth and D 2 W > in f2. Introduce local coordinates • • • , 0d-i) on 
S 1 ^ 1 satisfying the equality 



T w - . = (w xi ,. . . ,W Xd _ 1 ,-l) - (01,62,... ,9d-i,-\ 

y/1 + \VW\ 2 v ' \ \ 



d-1 

02 



Clearly, the first d—1 basis vectors of the ambient space constitute an orthogonal 
basis in the tangent space to S' d_1 at xq. Note that 

VW(x) 



where 9 = (6>i, . . . , #d-i), a = (xi, . . . , x d ) and 

d-1 

i=l 

Set W*(a;) = sup J/gf2 ((x, y) — W(y)). Since VW and VW* are reciprocal, one has 

x = VW*( — ===== ). 
Taking into account that Xd = W(x), one has 



h = (9, vw* (-_) ) - v/i-^F • ^(W* (-== 



Hence 

JT(0) = y/1- \9\ 2 W*( . 6 =) 

v ' v 11 V 1 - l^l 27 

on 7V(f2). This is equivalent to 

■v/TTR 



on Q. By approximation arguments these relations remain valid for every convex 
W in ft. Now assume that H is twice Alexandrov differentiable at 0. Clearly, 
H(0) = 0, V-ff(O) = 0. The same holds for W*. Using Alexandrov differentiability 
oi H, we get 

W*{x) = (^H{Q,-l)+ l -{{Dl)H^-l)x,x))+o{x 2 ). 

This means that D 2 a W* = H • Id + (Dg) a H. We get 1) by the duality relations for 
convex functions (see, for instance, The opposite implication follows by the 

same arguments. □ 

Clearly, if the surface is smooth and strictly convex, in the situation of the 
Lemma EOI one has 

K = -, — 5-^ . 

dct (H -Id + D 2 H) on 
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Definition 3.5. Let A be an arbitrary convex surface. We call the following 
quantity K "Gauss curvature of dA at x" 

(9) K{X): = det(H-ld+(D 2 g ) a H)on(x) 

if there exists a unique normal n(x), H is twice Alexandrov differentiable at n(x) 
and H ■ Id + (Dg)H is nondegenerate. 

If n(x) is unique, but H is not twice Alexandrov differentiable at n(x), we set 
K(x) — 0. The latter is equivalent to det D 2 a W{x\ 1 . . . , Xd-i) = if dA coincides 
locally with a graph of W : M d_1 -> E. 

Remark 3.6. Clearly, by Lemma \TM K(x) is well-defined for ^-almost all x, 
since W is % d_1 -a.e. differentiable. 



Remark 3.7. In the special coordinate system considered in the proof of Lemma 
one has K = dct D 2 W. Following the proof of Lemma |3"H1 one can easily understand 
that ((5J) holds almost everywhere in a non-smooth setting with F — W — Xd if the 
second derivative of W is understood in the Alexandrov sense. 

Recall an important result of McCann [25] , 

Theorem (McCann). (Change of variables formula for convex func- 
tions.) Let [i — f dx and v = g dx be two probability measures and V be a convex 
function such that v = // o VF -1 . Then for fi-almost all x one has 

g(VV)detD 2 a V = f, 

where D^V is the second Alexandrov derivative ofV. 

In the following proposition we deal with the Gauss map n: dA — > (non 
multivaled!) which is H d ~ 1 -a,.e. well defined. 

Proposition 3.8. For every A t — {x: <p(x) < t} the measure \ K ■'H d ~ 1 \dA t ^J °n _1 

is absolutely continuous with respect to H d ^ 1 and the following change of variables 
formula holds for every bounded Borel function f : S' rf_1 — » K: 

/ f(n)K dH*- 1 = f f dU d -\ 

JdA t Jn{dA t ) 

Proof. It is sufficient to prove this result for dA t n V instead of dA t , where V is 
a small neighborhood of a point xq G dA t with unique n{xo). Fix such a point 
and choose a coordinate system in such a way that n(xo) = ea and the surface dA t 
coincides (locally) with the graph of a convex function W : U C E rf_1 -> K, where 
U is an open ball containing and W attains its minimum at 0. In addition, we 
may assume that dW(\J) is a bounded set. Let U be a local chart of dA t PI V and 
parametrize a part of dA t f~l V in the following way 

U 3 (xi, . . . , x d -i) -t {xi, . . . , x d -i,W(x)). 

Since W is Lipschitz on U, the surface measure on dA t fl V can be computed 

in this chart by mo = (1 + IVW^ 2 ) 3 The Gauss map n is given by 

n = — 1 (-9 ai W,...,-9 ad ,W,1). 

This holds for almost every (x\, . . . , Xd-i)- 

li 



Identify the half-sphere S d 1 H{xd < 0} with its projection B d 1 on (xi, . . . x<j_i) 

and n with the mapping h: vw taking values in B^ 1 . Note that the 

li b ^i+|vw| 2 6 1 

surface measure on S^ 1 can be computed in the local chart 



(xi, . . . ,x d _i) (xi, . . . ,x d -i, y 1 -xf - ... - aj|_x) 

by mi = -O— U d -\ 

Note that n = F o WW is the composition of WW with the smooth mapping 

F(x) 
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which is nondegenerate everywhere. 

Writing the local chart expressions we get that the claim is equivalent to the 
equality = (K ■ m )\M + o h, where M + = {x 6 V^: detl^W 7 > 0}. 

Note that if is well-defined on M + . By Remark 13.71 we have 



K = det D a n = det D a W ■ det DF(x) o (VW). 
By the result of McCann the optimal transport VW pushes forward 

det D 2 a W-H d - 1 \ M+ 
to / H d ~ 1 \vw(M + )- Hence we obtain that the image of the measure 

K (1 + \VW\ 2 )iH d - 1 \ M+ = K ■ m \ M+ 
under VW coincides with 

det J F(a;)(l + |x| 2 )^ d - 1 |v M /(M + )- 
Now applying the standard change of variables formula we get that the image of 

\dW(M + ) 



detF(x)(l + \x\ 2 )^n d - 11 



under y = F(x) coincides with — ■ j ==='H d 1 |s(m + ) = m i\h(M + )- The proof is com- 
plete. □ 

The fact below follows easily from the McCann's theorem. 

Corollary 3.9. If V is a convex function satisfying det D a V\M — for some set 
M with X(M) > 0, then the image of X\m under W is singular to A. 

The proof of the following lemma can be found, for instance, in [7]. 

Lemma 3.10. If V is a convex function satisfying D 2 V > on M, then X\m ° 
W _1 is an absolutely continuous measure. 

We prove an analog of the McCann's theorem for the Gauss mass transport. 
We start with a change of variables formula for the mapping T -1 defined by a 
monotone-convex potential u. Since u and H are related by a smooth change of 
variables, it gives immediately a change of variables formula for H. 

Remark 3.11. We recall that u(z,r) (see Section 2) is convex in z and increasing 
in r. Hence one can define u r and (D 2 z ) a u H d -almost everywhere, where u r means 
a partial derivative of u in the classical sense. 
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Theorem 3.12. (Change of variables formula for u) The potential u satisfies 
the change of variables formula for T~L d -almost all (z, r) C R''" 1 x [0, R] 

u r ■ det (Dl) a u — 



Mr- 1 )' 

where 



Pi 



1 + z{ + . . . + z, 



i 

d-l 



d - x ' rzi,. . . ,rz d -i, 



Pi 



Proof. Fix an orthogonal coordinate system (xi , . . . ,x d ) and denote by e$ the corre- 
sponding basis. Recall that mapping T" 1 sends v = pxdx\^d x r 0tR -i to t L \r- 1 ({x d <o}) 
and admits a.e. the representation 

d-l 

T -1 = ^2 u Zi ■ ei + u*(V z u) ■ e d , 

i=l 

where u*(z) = sup^gRd-i ((x, z) — u[xy\ . Let us represent T~ l as the composition 
of two mappings T -1 = S 2 o Si, where Si : R^" 1 x [0,R] -> R^ 1 x[0,R] has the 
form 

Si(z,r) = (V z u,r) 
(all expressions are written in the Euclidean (z, r)-coordinates!) and 

d-l 

S 2 (z,r) = ^2zi- e~i+ u*(z,r) ■ e d . 

i=l 

Let us show that det \D^) u > almost everywhere. Indeed, set 

M: ={(z,r): Aet(D 2 z ) u = 0} . 

Assume that \(M) > 0. Then by Corollarv l3.9l and Fubini's theorem 

v = pi-H d - 1 \ M °Sr 1 

is a singular measure. Let us disintegrate v along the r-axis: 

v(r, z) = v z (dr) ■ fio(dz). 

Here no is the projection of v onto R d_1 and v z {dr) are the corresponding condi- 
tional measures. 

Denote by vq the projection of v onto (zi, . . . , z d —x)- It follows from the relation 
v o S2 = MlT _1 ({a; d <o}) that the image of Dq = (Ji/ Z (dr)) ■ (io(dz) under 

d-l 

(10) (zi,...,z rf -i) -> y^z t -Sj 

i=l 

coincides with the projection of t l \T~ 1 ({x d <o}) onto (xx, . . . ,x d -i)- Since the latter 
admits a Lebesgue density and (|10p is smooth and nondegenerate, one gets that 
Ho(dz) admits a Lebesgue density fo(z). Hence for 'H d_1 -almost all z the one- 
dimensional measure v z (dr)\s 1 (M) 1S singular. Note for T~L d ^ 1 -almost every fixed 
z the mapping r > u*(z,r) pushes forward v z (dr)\ Sl ( M j to a one-dimensional 
absolutely continuous measure. Since v z {dr)\ Sl ( M ) is a singular measure, one has 

u*(r, z) = —00 

(note that u* is decreasing in r) for T^ -1 -almost all z and ^(dr^nm -almost all 
r. This follows by duality from Corollary! 
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Next we note that 

M*(V z u,r) = — u r (z,r) 

A-a.e. Indeed, in the case of smooth functions with non-degenerated second deriv- 
ative it follows by differentiating the duality relation 

u*(V z u,r) + u= (V z u,z) 

in r. The general case easily follows by approximations. 

Thus u r = oo £>-almost surely on M. But this contradicts the assumption 
A(M) > 0. 

Since det(Z)5;) Q u > 0, by Lemma [3T0] the image v under S\ is absolutely contin- 
uous. Then one can apply the McCann's change of variables formula for 'H 1 -almost 
every fixed value of r. Applying the same theorem once again to S2 (for T-L d ~ 1 - 
almost every fixed z) one gets the result. □ 

Corollary 3.13. (Change of variables formula for H) Since T 1 and T 1 
are related by a smooth change of variables, one immediately gets 

H r - detjH- Id +(Dj) a H) 

Pi -Pol 1 ) —! 

H d -almost everywhere on Br. 

Corollary 3.14. The Gauss curvature K (x) = det(H Td+ [Dq)H) is well-defined 
and positive for ^.-almost all x. 

Recall that Dip denotes the generalized derivative of ip in the distributional sense. 
Since ip has convex sublevel sets, it is a BV function (see Q]). Hence Dip can be 
understood as a vector- valued measure satisfying 



/ (Dip, £) dx — — J ip div£ dx 



for every smooth compactly supported vector field £. We denote by \\Dip\\ the 
corresponding total variational measure and by |-D a <^| its absolutely continuous 
component and by |-D s y| its singular component. 

Theorem 3.15. (Change of variables formula for ip) The following change of 
variables formula holds for \x-almost all x G A: 

(KlD^-^ix) pi(T(x)) = po(x). 

Proof. Let A C A be a set, where K is well-defined and positive. By the previous 
corollary fJ,(A) = 1. Let us show that |Z> s (,o| (A) = 0. Indeed, otherwise we can find 
a set M s C A with X(M S ) = and \D s ip\(M s ) > 0. By the coarea formula for BV 
functions (see pQ, p. 159) 

< f K d\D s ip\ = f K d\\Dip\\ = [ [ K dU d ~ x dt. 

J M s JM S JO JdA t nM s 

By Proposition 13.81 and Fubini's theorem the latter equals 

f / t-^ dU d - l dt = f \y\-V-VdU d . 

JO JdB t nT(M s ) Jt(M s ) 

Since T pushes forward \i to v, one has A(T(M S )) = 0. We get a contradiction. 
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Applying again the coarea formula for BV functions we get 

f £(T)K \D a <p\dx = T f Z{T)K dH^dt, 

J A JO JdA t nA 

for any Borel bounded function £. By Proposition ^. 81 

/ ^T)KdH d - 1 = f m*H d -\y) 

JdA t nA JB t 

for almost all t £ [0,R]- Since T takes podx to p\dx, one gets 

m K\DM oT _ l Mdy= r,r dnd _ 1{y) ^ tl _ {d _ 1)dt ^ 

B r PO Jo V ifJ t / 

Hence, for /i-almost all y 6 B r , one has —iS^Xi T~ 1 (y) pi(y) = \y\~^ d ~^- The 
proof is complete. □ 

Corollary 3.16. Comparing different change of variables formulae, one gets 
p-almost everywhere. 

4. SOBOLEV ESTIMATES FOR (f 

The main goal of this section is to establish some natural Sobolev estimates for 
if (Theorem I4.5|) . The proof is based on the integration- by-parts and change of 
variables formulae. 

Before proving Theorem 14. 5 1 we establish some |Vc/?|l°o -bounds with the help of 
the classical maximum principle. These estimates have an interest in their own, 
they will also serve as an intermediate step in Theorem 17.11 

It will be assumed below that H r , Hg , Hgg are continuous and continuously dif- 
ferentiable in r (except, maybe, the origin) up to the boundary. We also assume 
without loss of generality that H > and H(0) — (this can be achieved just 
by shifting A and assuming that ip(0) = 0). The estimates obtained below do not 
depend, however, on higher derivatives (see in this respect Remark I4.4[) . 

Let us set 

P = Pl r d - 1 . 

Since H is smooth, it satisfies 

P = po{T- v ) H r ■ det(# • Id + D 2 e H) 
up to dB R . We recall that H r = l/|V( y 3(T- 1 )|. 
Proposition 4.1. a) Assume that for some C > 

|Vp |<Cpo + ^, P<C, 
and there exists u: (0, R] — >■ R with u € L 1 ([a, i?]) for every R > a > such that 

—— < u\r). 
p - V ) 

In addition, assume that dA is smooth, Xq = inf^g^ K{x) > and PqIqa < C, 
P\dB R > Then 

r R 



H r >Diexp(—J u(s)ds). 



is 



In particular 



\V(p\ < £>2exp( / u(s)ds 



with Di, D 2 depending on d, C, Ao, R- 
b) Assume that for some C > 

|Vp | 



and 



<C, Po >-, P<C 
Po C 



In addition, assume that dA is smooth, A = sup xedA K(x) < oo, po\dA > C, 
< A. 

R—C 



P\dB a < h and 



H> er 



for some e > 0. Then 



with Di,D 2 depending on d, C, R, e, A . 

Proof, a) We are looking for the minimum of H r e? on Br \ B ro r > 0, where 
/ = f{r) is a radially symmetric function to be chosen later. Assume that the 
minimum is attained at some point x OBr . We deal with the local coordinate 
system (r, 9) as described at Section 2. Let us differentiate \ogH r + f(r) along r 
and every 6i at this point. One has 



0, ^ > -fr- 

£l r 



The second differentiation yields H r g i9i > 0. Rotating the coordinate system when 
necessary we may assume that Dg H is diagonal at x - Differentiating the change 
of variables formula in r yields 



P r _ (VpoiT- 1 ), H r -n + ££1 H r9i ■ (*) , H, 



Hence 



> Hr 

> H r 

This implies 



Po(T-i) 
(WpoiT- 1 )^) 



+ 



+E 

i=l 



(VpoiT- 1 )^) 



fr + H r 



? = 1 



H + H, 



OiBi 



(VpoiT- 1 )^) 



-f r + (d-l)H r [ 

- f r + (d- l)H? 



det(H + D 2 efj H) 



-J- i 

W 1 < 



d-l 



P 



Po(T-i) 



Pr 



+ .fr 



H r 



(V P0 (T-'),n)p-^(T-i) 
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Applying Holder's inequality one gets 

d 

H*- 1 < 



1 


f P ~ 


1 

d-1 




d-l 









2 



P^\V ■ po (T~ 1 )\p d - 1 (T- 1 ) 



where C\ depends only on d. Let / be of the type 

/ = —C%r — I u(s)ds. 



One gets 
1 



-A- Co 

-H*- 1 < =2- 

2 



P 



-Po(T^) 

Then it follows from the assumption of the proposition that the right-hand is neg- 
ative for a sufficiently large C2 > 0. This contradicts the estimate H r > 0. 
This means that 



H r exp(— CyX 



u(s)ds) 



can attain its minimum only at OBr- Taking into account that 

H r \dB R >C 3 {C,R) inf K{x) 

X^oA 

one gets the desired estimate. 

b) In the proof we use an idea from [34] . We are looking for the maximum of 



H - g{r) 

on Br \ B ro , where g = |r. Note that H — g > |r. Assume that log H r — \og(H — 
g(r)) attains its maximum at xq with |xo| < R (otherwise the estimate is trivial). 
Then at this point 

He, 



H rr H r — g' ^ H r g 
H r H — g H r 



H-g 



0. 



The second differentiation gives 



H r ~ H - g 

Differentiating the change of variables formula one obtains 
P H r f-^ H + Hg iSi p (T l ) 



Hence 
d-l 



d-l 



E H r + H r g i g i H r 
H + He." ~ H -n 2^ 



H-g + Hg t9i H, 



d-l 



<(d-l) 



H-gfr! H + Hg 



= ^-g{ d - l -3Y [ WTH-g 



Hr 



H-g 
H, 



1-5 



(d-l)^(l-, 



7/ r 



Po(T- 
P 
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Next using 

d-l 



»=i 

we get 



^ ((VpoiT- 1 ), n) + ^J— ^ if 9i (VpolT- 1 ), ei )) . 



i=i 

Taking into account the assumptions, boundedness of Hg i and H, we get 

-^(VMr->),r-.)< Cl ^. 

Thus we obtain 



P H — g H-g\ * V P > H-g 

Multiplying this inequality by H — g, using the assumptions of the theorem and 
boundedness of H we get 



-C(H -g)< C 2 H r -g'-(d- 1 :■/,//, d ■ r "' ' ' ' 



< C,^-^- + ^ f- 5 ' + C 2 H r ) < 9l±^L < W-D + c 7 (-)' 
g q \ J r 2 \rJ 



Thus implies 

{H — g) | d ( j i ^ u \ ^ C5 + C§H r 
9 

Hence 

\H-g) ~ r d+n^' 

This gives the desired result. □ 

Remark 4.2. The proof of a) can be generalized to the case of pre-limiting potentials 
H t (see Section 2). Since the computations are quite involved, we give only some 
intermediate results. For simplicity let us skip the index t and write H instead of 
H t . Choose a function / in such a way that f(x) <~ — (d — 1) ln(r — r ) + for a; close 
to dB ro and assume that the minimum point x does not belong to OBr. One has 

d-l 

T- 1 = (H+ j^H r ) ■n + Y.He,- 



The derivatives of T 1 at x satisfy 

i + 2 r d ~ 1 

T" 1 = (—H r + j^jH rr ^ ■n + Y.Hrer ei. 

Choosing an appropriate basis, we may assume without loss of generality that 



H 0i0j = 

18 



d-1 



for i 7^ j . Then 

z— 1 

t^n([iTi a+ f^] + M- 



+ 



i=l i^j 
At the minimum point one has 



(11) 
(12) 



11 r 



+ (/" - (ff) > 0, H r e i6i > 0. 



The reasoning from the above proposition leads to the following estimate: 

(Pi)r | >H ft + 2 r ^(n.VpoCT" 1 )) 

pi r 



+ 



r U + l t+1 1 ) poiT- 1 ) 
(t + 3)/' + r((/') 2 -/") 



(t + 2) - rf 
tit 1 (, rf 



y ' r \ t+l)Wt+l 1 + tJ 



r \ t+1, 
Choosing an appropriate / one gets the desired bound. 

Corollary 4.3. Assume that 

P<C,^< Po , 
dA is smooth and uniformly convex, 

P 

< C 



1 

d-1 


-poiT- 1 )- 


i 




- Pi - 





Vp 




Pr 


Po 




P 



and p \oA < C, ± < P\ dBn . Then Dnp d < |V^| < D 2 for some Di, D 2 > 
depending only on d,C and dA. 

Remark 4.4. We have proved the above estimates assuming smoothness of H. But 
the hnal results do not depend on the bounds of the derivatives of H. We give 
some sufficient conditions for H to be smooth in Section 7. Applying smooth 
approximations it is possible to show that the estimates remain true without extra 
smoothness assumption of the solution. In particular, the upper bound on |V</j| 
implies the absence of a singular part for Dip. 

Theorem 4.5. Assume that pi = -§=r- Then for every p > there exist C p .r > 
such that 



(13) c PtR f iv^r 1 dp < f 

J A J A 

(14) C P , R f |V^| P+1 dp < f 

J A J A 



L 


Vp 


p+i 




Pa 




f 




P+i 


A 


Po 
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dp,+ f \V<p\Ppo dH d -\ 

JdA 

dp+ f K- p p p Q +1 dU d -\ 

JdA 



Proof. Under assumptions of the theorem the change of variables formula reads as 

CK\Vip\ = p . 

Computing DT is the standard frame {n, vi, • • • , v<i-i}, we get 

OT -(I*|H » 

where 

The Jacobian matrix of S = T^ 1 computed (r, 6) coordinates has the form (recall 
that d n — d r and = rd Vi ) 



DS = 

with 



H r 
c* B 



c=(H rei ), B=[ — 

Recall that H r (T) = r^t- Since DS(T) — DT' 1 , one can also assume that A and 
B are diagonal (at a fixed point). Denote by Xi the eigenvalues of A. Then using 
DT o DS(T) = Id one easily obtains 

^ + H,,<T)X,=0. 

Next we find 

Vnn = 0n|V<^| = d n (l/H r (T)) = -^^(^H rr (T)\V<p\ +J2 H re i {T)(d n n,v i ) 

1 



i_(iI rr (T)|V^-X: 



1 ,„2 



^rWj \ . H rr (T) 



Taking into account that tp has convex level sets (hence div-r^r > 0), we get 



\V<p 

div(^-^|V^) > |V^ +1 +p V \V V \ p -\nn > |V^ +1 



|V<p| / iJ 



o 1 . 



Thus 

(15) |V^ +1 <div(^|V^)+^ 
Integrate (p~5|) over v4 with respect to /i. One obtains 

; d iv(^|V^)po d X = R |V^|> - J a V^^m* *> 

Vp0 P+1 d». 



oT. 



<R / \V<p\ p po dU d - 1 +e / \Vip\ p+1 d^ + N(e,p) / 

J OA J A 



P0 



2(1 



Applying the change of variables and integrating by parts we get 



P 



f rH rr f rH rr p f 

T7fn oTd ^^P 77^2 d " = '77Tl (VH^*~ ,x)pi dx 

JAtir JB R rlr P + 1 J B R 

= — ^r/ H-*- l p 1 dU d - 1 + j-?— [ H-^ Pl dx 

P + 1 JdB R (P + 1) Jb r 

f |V^| P+1 dn--^r ( |V^| p p dU d -\ 



OA 



(P + 1)Ja P+ 1 
The obtained estimates imply immediately (|13p . Estimate (|14p follows from (fT5|) 
and the change of variables formula. □ 



Remark 4.6. Estimates of these type are also available for the pre-limiting poten- 

/i \ 

oT 



tials. For instance, for T = ip-S^r\\jLp\ w one has 



for p — 2 + — ^ . Then one can show that for q > 

/ ¥>|VvlVo dH 1 + / — ?+ V +1 d/i > C g / |V(^|« +1 ^. 

J dA J A PO J A 

Remark 4.7. The result can be easily generalized to the general case of a continuous 
rotational invariant density v = p v dx = p v {r) dx. 

Indeed, take a mapping T sending v to — and having the form T(x) = f(r) — . 
The function / satisfies 

rp v {r) = f'(r). 

Note that 4>]^!f\ ; where ip = fif) sends pb to Applying (fT3| to ?/> we get 
C / ^ +1 p£ +1 (^)|V^ +1 d/x < / ^ P+1 dp+ [ \Vif\ p+1 po du d - 1 . 



OA 



Remark 4.8. It looks possible to prove L°°-bounds on \V<p\ using the parabolic 
maximum principle (see the next Section) and assuming high integrability of |Vpo|- 
Estimates of this type for the potential u have been obtained in [20] . Results from 
|20) are not directly applicable to our situation, since we need to consider u in 
unbounded domains. 



5. Variants of the parabolic maximum principle 

For every convex V we denote by \dV\ (B) the associated Monge- Ampere measure 
of the set B, which is defined as follows: 

\dV\(B)=X({\J dV(x)}), 

x€B 

where dV is the subdifferential of V at x. 
For smooth V one has 

dV = detD 2 V dx. 

This means that W sends dV to Lebesgue measure if det D 2 V ^ 0. 

Recall that for every continuous function / on a convex set A one can define its 
convex envelope /* which is the supremum of all affine functions less than /. The 
set Cf — {x: f(x) = f*(x)} is called the set of contact points of /. 
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According to the elliptic maximum principle (also called Alexandrov maximum 
principle or Alexandrov-Bakelman-Pucci principle) every continuous function / on 
a convex set A c K d satisfies 

sup / < sup / + C ■ diam(A) [#/*(£/) 

A OA L 

whcre C depends only on d. If / is twice continuously differcntiable, this implies 
sup / < sup / + C ■ diam(A) [ / | det D 2 f{x)\dx\ d , 

A OA k/.D 2 /(x)<0 J 

where C depends only on d. Equivalently, passing to g — sup A f — /, one gets that 
for every non-negative g 



inf q < C ■ diamM) 

dA* - X ' 



det D 2 g(x)dx . 

lD 2 g{x)>Q J 

A parabolic version of the maximum principle was obtained by Krylov (see [25]). 
Later Tso [35] simplified the proof in some special cases and gave extensions in 
some particular cases. 

In this section we prove some other variants of the parabolic maximum principle. 

Definition 5.1. For a continuous function / defined on a convex set A consider 
its sublevel set A t = {/ < t} and the convex envelope conv(A t ) of A t . Every point 
x e Int^4 satisfying x £ A t l~l <9conv(A t ) for some t we call a contact point of A t . 
The set of all such points will be denoted by Cfj. 

We denote by Sir 1 the upper half of the unit sphere in R d . For every set 
f2 = |(r, 6): R\ < r < i?2, 6 £ Qj, where Q C Si^T 1 is a spherically convex set, 
we denote by 

dptt = Q x R 2 U dQ x [R U R 2 ] 

its parabolic boundary. 

Theorem 5.2. 1) Let v be a twice continuously differentiable function on a 

convex set A C K d . Then there exists a constant C = C(d) depending only 
on d such that 



(16) 



supw(a;) < sup v(x) + C(d) 

iEi x£dA 



| | if da;. 



where K(x) is the Gauss curvature of the set d convjy: v(x) < v(y)} at x 
2) Let be a set of the type 



{(. 



Ri<r<R 2 



0} 



(17) 



with a spherically convex Q C satisfying dist(Q, dS^ 1 ) > 0. Then 

for every twice continuously differentiable function f : Q — > K satisfying 
suP^eSpO / > 0, one has 

\f r det(f -Id + D 2 J)\ 



sup / < Ci sup / + C 2 



-dx 



r. 



where T f = {x £ ft : f r < 0, / • Id + Dgf < 0}, and constants C\,Ci > 
depend only on d and Q. 
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Proof. 1) Set / = (M - vf/ d , where M = sup^ v. The estimate (fTO)) is equivalent 
to 

(18) mif d <cf f d ~ x \Vf\Kdx. 

9A J Cfl 

For every < t < inf qa f let us consider the set A t = {x: f(x) < t} C A and 
its convex envelope conv(A t ). Since A is convex, conv(A t ) lies inside of A and, in 
addition, dist(conv( J 4 t ), dA) > 0. Set: Ct — <9conv(^4 t ) n A t . Since / is smooth, 
the image of Ct under the Gauss map n of 9conv(^4t) covers the unit sphere. Hence 
the image of 

U Ct=C } , 

0<t<inf 8 A / 

under T = / ■ n coincides with {x: \\x\\ < mi dA f}. One has det DT = f d - 1 \Vf\K. 
The result follows from the change of variables formula. 
2) Let us consider the set of vectors V satisfying 

a) (v, n) < M. for all x G O 

and 

b) (v, n) > m, for all x G <9 p £l, 
with n = A, M = sup xeA /, m = sup^g^ /. Since dist(Q, dS^ 1 ) > 0, the set of 
vectors v satisfying b) is non-empty and has the form 

(r,9): r > C(Q)m, 6eQ 

for some set Q C S^ 1 and a constant C(Q) depending on Q. If M < C(Q)m, the 
claim is proved. If not, then V is nonempty. Consider the set 

B = {(r, 9) : C{Q)m <r<M, 0eQ}cV. 

Clearly, 

C (Q)(M - C(Q)m) d < X(B). 

It remains to estimate X(B). For every a € B define M a — {x: f(x) — (a,n)}. 
Conditions a) and b) imply that M a is non-empty and contained inside of f2. Hence, 
there exists a point xq € M a in the interiour of J7, where |x| attains its maximum. 
One has at this point 

f(x ) = (a, n), 

fv(x ) = (a,n v ) = - — Aa,v). 

\x \ 

for every unit w_Ln. This implies that 

Def = a - (a, n)n. 

In addition, 

f r (x ) < 0, D 2 9 f(x ) < D 2 e (a, n) = -(a, n) • Id = -f(x ) ■ Id. 
Hence BcT/. Set: 

d-l 

S = f(x)n + \x\J2 /vi (x)vi = f(x)n + D g f(x). 

i-l 

Note that S(xq) — a. This means that 

S(T f ) = B. 

23 



By the change of variables formula 



X(B)<I d etDSdx=l \fr**{f-U + Djf)\^ 



The proof is complete. □ 

Remark 5.3. Inequality (jTTJ) implies a form of the parabolic maximum principle 
(see [35]). Assume that sup 9pn , / = 0. Then 

|/ r det(/-Id + Dg/)| 



(19) sup/<C 



-dx 



Set u = vl-|- z 2 f, where z and x are related by the change of variables described 
is Section 2. Using 

det(/ ■ Id + Dgf) = det((l + z 2 f/ 2 D 2 z u^j 
and trivial uniform estimates one gets 

(20) sup u < C{d, Q)\ \u t - det D 2 x u\ dtdx 

r« = {ut < 0; D 2 u < 0}, for any u with sup ap Q u = and a cylinder SI = [0, R] x Q 
with convex Q. To remove the restriction supg^ u — one applies the estimate to 
u = v- sup dpil v. 

Remark 5.4. The above variants of the parabolic maximum principle are naturally 
related with transport mappings of the type 

T = cp^L, S = H-n + D e H. 

Both variants of mappings can be obtained from the "elliptic" transportation 
by scaling procedures (see Section 2). The transportation by gradients arc natu- 
rallly associated with the elliptic maximum principle. Is it possible to derive both 
parabolic maximum principles from the elliptic one? 
1) Elliptic maximum principle implies (|16p. 

We prove that for every continuous / > on a convex set A C M. d satisfying 
mi xe A f(x) = and every < p < 1 there exists a constant C = C{d) depending 
only on d such that 

inf < C diam*(A) \dW*\ o S^fcw), 

where Cw is the set of contact points of W = and S p (x) — \ x \ i- P ■ In 

particular, if / is twice continuously differentiable, one has 

(21) inf/ d < 1+p ) < C diam dp (A) / det Dff .J 7 * )dx. 

Clearly, letting p — > we deduce an equivalent form of P^)l from (|2T|) . 

Proof: Let xo be a point satisfying f(xo) = 0. If xo € dA there is nothing to 
prove. Thus we assume that xo (fc dA. Let V be the convex function whose graph 
is the upside-down cone with vertex (xq, 0) and base A with V = m on dM, where 
m = inf Ke si ^if 1+ ^( x )- It is easy to check that 

B m/ <tom [a) CdV(x )cdW.(C w ). 
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Note that the measure with density 

p: =det£> 



is the image of Lebesgue measure under S p . Hence 

J pdx < IdW*] o S- 1 (c w y 

s ,n/diam(f!) 

The direct computation yields p = p r d ( p ~ x \ This immediately gives 

m/diamcn) 

with C depending only on d. This proves the first part. 

Finally, (|2"Tj) can be obtained by direct computations. We just notice that 

{x: D 2 f 1+ ^{x) > 0} C C w - 

2) Elliptic maximum principle implies (|1T[) with SI = Br and symmetric 

/• 

Let /: Br -^Ibea symmetric (/(— x) = f(x)) bounded function. Assume that 
/ is twice continuously differentiable at every x ^ and mi x£ B r f( x ) < 0. Then 
there exists a constant C = C(d) depending only on d such that 

/ r det(/-Id + Z?g/) 



(22) inf / < C 

oiSu 



dx 

r-/ r 



Proof: For every t > consider 

W t {x) = \x\ /(x|xpTTT) 

defined on B R i+t . One has 

R 1+t inf / = inf w t (\z\ 1+t ) = inf w t . 

dB R zGdB R dB Rl+t 

Since u>t(0) = 0, by the elliptic maximum principle 
'inff. 



Indeed, Wt is twice continuously differentiable everywhere in B r except, maybe, the 
point x — 0. Without loss of generality one can assume that inf b r w t < 0. Since 
Wt is continious, inf w t is attained at some point x. Since Wt is symmetric, the 
points (x,w t (x)) and (— x,w t {x)) belong to a horisontal supporting hyperplane to 
the graph of wt- Since u>t(0) = 0, clearly ^ C Wt . This justifies the above estimate. 
Set: S t (y) =y\y\ t . Then 

/ detD 2 w t dx = det D 2 w t (S t ) det DS t dy 

det D(Vw t {S t ))dy. 
Direct computations yield 

rf r 
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Hence 



limdet£'(Vw; t (S , t )) 



f r det(f + D 2 g f) 



„d-l 



The proof is complete. 

3) Does elliptic maximum principle imply (|20p ? There are good reasons 
to believe that the elliptic maximum principle implies (|20[) . This problem seems 
to be rather involved technically and we do not consider it here. We just give a 
proof in a particular simple case. Let / satisfy all the assumptions from item 2). 
In addition, assume that / = C outside of ft' U (—ft'), where 

ft' = |(r, 6») : < r < R, 6 e Q'j 

and Q' C S^ 1 satisfies dist(Q', dS^' 1 ) > 0. By the previous result 

r frdet(f -Id + D 2 e f) 

Arguing as in Remark 15.31 we get that 



inf / < C 

aw - 



„d-l 



-dx 



inf u < C(d, Q) 



\ut ■ det D 2 x u\ dtdx 
satisfying 



holds for any bounded u: (0, R] x R d_1 , 

1) u is smooth on (e,R) x R d_1 

2) u is constant outside of ft — (0, R] x Q with convex Q C 

3) info u = 0. 

Passing to u = sup^ ip — ip we obtain 



pd-l 



sup <p < C(d, Q) 
n 



\<pt ■ det D 2 x Lp\ dtdx 



for any smooth compactly supported ip with supp(</?) C ft. 

6. ISOPERIMETRIC INEQUALITY 

We discuss two apparently different proofs of the isoperimetric Euclidean inequal- 
ity for convex sets (it is well-known that the general case can be easily reduced to 
the convex one). First of them due Gromov. It is worth mentioning (this was 
pointed out to the author by S. Bobkov) that arguments of such type go back to 
Knothe [24]. More precisely it has been shown in [24 that the Brunn-Minkowsky 
inequality can be proved by transportation arguments with the help of triangular 
mappings. The second proof comes from the differential geometry. Our aim is to 
reveal a remarkable similarity between probabilistic and geometrical points of view. 

1) (Mass transportation. Probabilistic approach.) We follow the mass 
transportation arguments but use the Gauss mass transport instead of optimal (or 
triangular) one. Let A C M. d be a convex set and T = ^j^y send Lebesgue measure 
on A into Lebesgue measure on Br, where Br is a ball of the same volume. By 
the change of variables (see the previous section) 

<p d - 1 \D a <p\K = l. 

Hence by the arithmetic-geometric inequality 
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1 = det D a T < 



where D a T is the absolutely continuous part of the distributional derivative DT. 
Clearly, 

X(A) < / div T dx — -JL / (T,n> dH^ 1 < —^—% d ^ 1 {dA). 

d-1 J A d-1 J dA d-1 

Taking into account that X(A) = X(Bji) = CdR d , one easily recovers the classical 
isoperimetric inequality. 

2) (Curvature flows. Geometric approach.) The same proof can be rewrit- 
ten in the language of curvature flows. The curvature flow proofs are well-known 
in differential geometry (see partial results on the Cartan-Hadamard conjecture in 
[33] . [30]). Let A t = {x: (f(x) < t}. For convenience we assume that <p is smooth 
on {x: ip(x) > 0} (which is indeed the case for smooth strictly convex dA). Note 
that A t are expanding with the speed y^r^- The enclosed volume X{A t ) evolves 
with the speed which can be exactly computed by the Gauss-Bonnet theorem 

^X(A t ) = t d - x [ K dU d - x = t^H^iS 11 - 1 ). 
at JdA t 

Hence 

X(A t ) = X(B t ). 
In the other hand, it is known that 



— U d - X {dA t ) = t^ 1 [ KH dU a 
at JdA, 



I dA t 

where H is the mean curvature. By the arithmetic-geometric inequality K x l^ d ~~ x > < 
-^j. Hence 

^H^idAt) >(d- l)^" 1 / dH*- 1 . 
at JdAt 

By Holder's inequality 

i/(d-i) 
i/(d-i) 



>(d- i)t d - 1 ( k dw 1 - 1 ) 3=1 (n d -\A t yj 



= (d - i)t d - 1 (u*- 1 ^- 1 )) 3=1 (u d -\A t )) 

Integrating in t one obtains 

n d ~ l {A t ) > t d n d - 1 (s d - 1 ) = u d - x {B t ). 

The proof is complete. 

7. On Holder's regularity of the Gauss mass transport 

The elliptic and parabolic Monge-Ampere equations belong to the family of the 
so-called fully nonlinear PDE's. See [12] (and [19] for the special case of the Monge- 
Ampere equation) . A short survey |26j presents the developments of the main ideas 
of the nonlinear PDE's theory. 

The connection between the variational Monge-Kantorovich problem and the 
elliptic Monge-Ampere equation was revealed by Brenier (see [36]). In [34] the 
existence of the Gauss curvature flow for smooth data was established by solving 
the corresponding equation of the parabolic Monge-Ampere type. 
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Contributions to the regularity theory of the elliptic Monge-Ampere equation 
were made by many authors, including Alexandrov, Calabi, Yau, Pogorelov, Krylov, 
Spruck, Caffarelli, Nirenberg, and Urbas. There are several approaches to the reg- 
ularity theory of nonlinear equations. A classical one is based on differentiating of 
the underlying equation. Taking the second derivative one obtains another equa- 
tion which is linear with respect to higher derivatives. Then one applies a priori 
estimates from the linear theory. This was a common way for studying the nonlin- 
ear PDE's before the results of Krylov, Safonov, and Evans on a priori estimates 
for nonlinear uniformly elliptic operators. See [25] for details. 

Unfortunately, the elliptic Monge-Ampere operator 

u — > det D 2 u 

is not uniformly elliptic even in the class of convex functions. This is the reason 
why the Krylov-Safonov-Evans theory is not applicable directly. The regularity 
problem for the elliptic Monge-Ampere equation was solved in sufficient generality 
by Caffarelli. Combining the nonlinear regularity theory and deep geometric consid- 
erations he proved, in particular, that the solution V of the optimal transportation 
problem 

g(VV) dctD 2 V = f 

for probability measures / dx and g dy with compact supports X and Y is (2 + a)- 
Holder continuous inside of X provided /, g are Holder continuous, bounded away 
from zero and Y is convex. 

Many regularity results for the Gauss curvature flows (see [34], [3]) were obtained 
by using the classical way of differentiating the evolution equation. Similar to the 
the elliptic case, the parabolic maximum principle plays a crucial role in the study 
of this problem. 

A parabolic analog of regularity theory for uniformly operators has been devel- 
oped in [40] . 

The regularity of the parabolic Monge-Ampere equation was studied by Krylov 
[25] . Ivochkina, Ladvzhenskava [23] . Gutierrez, Huang [20], R.H. Wang and G.L. Wang 
[38] . [39] (see [26] for references). Some interesting results were proved by prob- 
abilistic methods (optimal control and stochastic differential equations), see [2B], 
[32] . A parabolic analog of the Caffarelli theory for the elliptic Monge-Ampere was 
developed by R.H. Wang and G.L. Wang in [38], [39]. They studied the parabolic 
Monge-Ampere equation 

(23) u t det D 2 z u = f(t,z) 

on the domain Q = O x [0, T] with given values u = ip(t, z) on the parabolic 
boundary d p Q. It was shown in 38 that under the assumptions that 

1) is compact, strictly convex with C 2 -boundary 

2) / is positive Lipschitz continuous on Q 

3) if € C 2,1 (Q) with ip t > 0, D 2 z Lp > on Q 

there exists a solution u G C'fJ c a ^ 2 ' 2+a for some a > 0. A measure-theoretic 
interpretation (the parabolic Monge-Ampere measure) was given in [39) . 

For further generalizations and refinements, see [28 . Sobolev estimates for ([23]) 
are obtained in [2Tj] . 

We analyze below the regularity result of Tso. The reasoning from [M] can be 
easily generalized to our situation with the help of our results from Section 4. We 
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will not repeat the lengthy reasoning from [33] and give just a brief sketch of the 
proof. 

Let C k,a (Bft) be the parabolic Holder norm (see [34]) on functions 

f(r,d): [0,R] x S^ 1 h> R. 

It was established in [33] that for p = 1, p x = ^(0) G C 2+Q (S' ci " 1 ) with some 
a > and every i? > r > there exists C such that 

|-ffrlc/3/2,3(_B R \ Bro ) + \DeH\cf>/ 2 >P(B B \B T0 ) < ^ 

for some /3 > 0. 

Using the estimate < c ra < £T r < C ro for smooth ff from Section 4 and 
repeating the arguments from [34] it is not hard to verify Theorem 17.11 below, 
which is a generalization of Theorem 4.2 from [33] • Clearly, a solution H obtained 
in this theorem coincides with the potential H of the corresponding Gauss mass 
transport by the uniqueness theorem from [9]. 

Theorem 7.1. Assume that p\ G C 2 ' a (Bp), p G C 2,a (A), A is uniformly convex 
and Ha G C 2,a (S' d_1 ). XTien a solution H to (OJ) wif/i i?|a^4 = i?4 exists. In 
addition, 

HeC 4 > e (B R \B ro ) 

for every ro and 

( 24 ) l^lc' 3 / 2 " 9 (S n \B ro ) + \DoH\cl>/2,t>(B R \Br ) < ^ 

holds for some positive /3,e, C depending on ro, R, the curvature of dA, the Holder 
and uniform bounds on po, pi- 

Sketch of the proof: One proves the existence of a solution to ©. 

1) The classical short-time existence result implies that a smooth (C 4 ' £ ) solution 
to ([6]) with a given initial value H(R, 9) = Ha(0) exists for t G [R — e, R] (see, for 
instance [16], Theorems 2.5.7, 2.5.9). Let [R*, R) be the maximal existence interval. 
Assume that R* > 0. Applying the change of variables formula, let us estimate the 
volume enclosed by the hypersurface determined by H(R* , 9). One concludes that 
there exists a sphere contained in all hypersurfaces determined by H(r,0), r > R*. 
Taking the center of this sphere as the new origin one can assume without loss of 
generality that H is strictly positive on [R* , R] . 

2) The results of Section 4 give ^ < \H r \ < C for some C > and every r G 
[R*, R], Following [34] one obtains that A > H + D 2 H > A for some constants < 
A < A, t € [R*, R}. This can be shown by differentiating twice the equation in 9 and 
applying the classical maximum principle to a suitable function (see also Pogorelov- 
type arguments in 18!, Theorem 17.19). Alternatively, one can use Caffarelli's 
result [13 on bounds for principal curvatures of a smooth convex set under the 
assumption that the corresponding Gauss curvature is positive and bounded. 

3) Differentiate §6§ in r. The Krylov-Safonov estimates (see [25]) imply that the 
parabolic Holder norm of H r on [R* , R] is under control. The same holds for Hg . 

4) It remains to prove Holder's continuity of H + DgH. The arguments follow 
|34j . Let us indicate the main difference. To estimate oscillation of Hgg (or u zz ) we 
need an estimate for the additional term 

| logpoCT" 1 )^) - logpoCT- 1 )^)! < sup ^IT- 1 ^) - T~ 1 (y)\. 

Po 
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(See [31], Theorem 4.1, (4.3)- (4.4)). Then we estimate T~ 1 (x) - T~ 1 (y) by a 
parabolic Holder norm of Hg (see item 3)). Thus we get 



H<=,i 



(^1 + 5/2,2 + 6 



< c 



1 



for some 8,e > 0. Then the parabolic interpolation inequalities (see [27], Theorem 
8.8.1.) complete the proof. 

5) Since we have managed to keep control on the norms of derivatives of H on 
[R* , R] , the solution exists for r < R* by the short-time existence theorem. Hence 
R* = 0. The proof is complete. 

This work was supported by the RFBR projects 07-01-00536 and 08-01-90431- 
Ukr, RF President Grant MD- 764.2008.1, and the SFB701 at the University of 
Bielefeld. 
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